Non-perturbative saddle point for the effective action of disordered and interacting 

electrons in 2D 
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We find a non-perturbative saddle-point solution for the non-linear sigma model proposed by 
Finkelstein for interacting and disordered electronic systems. Spin rotation symmetry, present in 
the original saddle point solution, is spontaneously broken at one-loop, as in the Coleman- Weinberg 
mechanism. The new solution is singular in both the disorder and triplet interaction strengths, and 
it also explicitly demonstrates that a non-trivial ferromagnetic state appears in a theory where the 
disorder average is carried out from the outset. 
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Understanding the combined effects of disorder and in- 
teractions in electronic systems has proven to be both an 
extremely interesting and difficult problem. The issue 
has regained a new surge of interest since the discovery 
by Kravchenko and co-workers Q of a possibly conduct- 
ing state in two-dimensional (2-D) Si-MOSFETs (sec also 
Ref . 0] ) . Although there has been intense debate within 
the theoretical community as on the origin of the tran- 
sition H , little has been accomplished that is in as solid 
grounds as the scaling theory of localization for the non- 
interacting problem [Q. 

The initial attempt to establish a scaling theory for 
the interacting problem was put forward by Finkelstein 
B, who studied an extension of Wegner's non-linear 
sigma model containing singlet and triplet interaction 
couplings. The renormalization group (RG) flow takes 
the interaction coupling constants to strong coupling, 
away from the perturbative starting point of a diffu- 
sive Fermi liquid state; hence no conclusive picture has 
emerged from this approach. Among the outstanding 
theoretical issues is the nature of the magnetic state sig- 
naled by the divergence of the triplet interaction. 

Let us begin by briefly discussing the breakdown of 
the RG flow for the interaction. The one loop RG equa- 
tion for 72, the ratio between the triplet and the singlet 
interactions, is given by 
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where t is the resistance. Although the resistance also 
gets renormalized, it is instructive to solve this equation 
for a "fixed" i, which leads to a breakdown at an RG 



the flow equation for the BCS interaction V < 0: 

dl 
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This equation also breaks down at a finite length scale 
- the coherence length. At this length scale, the non- 
perturbative physics of pairing takes over. Similarly, 
the breakdown of Eq. ([l]) implies that non-perturbative 
physics dominates the behavior of the system. Such 
physics cannot be accessed perturbatively from the dif- 
fusive saddle point. 

However, the main difference between these two ex- 
amples is that, in the case of BCS superconductivity, 
we know what the physics of the non-perturbative fixed 
point is from the BCS mean-field solution. Shankar's 
RG procedure allows us to find that there is an insta- 
bility of the Fermi liquid, but it alone does not access 
the non-pertubative BCS solution. The situation is com- 
pletely analogous in the diffusive Fermi liquid problem: 
we know that the non-interacting saddle point is unsta- 
ble under the RG flow, but we cannot determine the non- 
perturbative physics solely from the flow. 

In this paper, we find a non-perturbative self- 
consistent solution for the saddle point of the interacting 
and disordered electronic problem by looking at the one- 
loop effective potential of the non-linear sigma model. 
The new saddle point spontaneously breaks spin rotation 
symmetry for any value of the triplet coupling constant. 

The starting point of our calculation is Finkelstein's Q- 
matrix model of unitary class Q . The disorder-averaged 
7V-replica partition function of the interacting problem 
reads 



scale I = t (i+ l2 ^ , or at a length scale A = A 1 e'C+^i (A 
is a momentum cutoff). 

It is very useful to draw an analogy to the problem of 
BCS superconductivity at this point. Within Shankar's where 
renormalization group approach to fermions Q , one can 
perturb around the Fermi liquid fixed point, and obtain S[Q] 
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with the tensor 7 containing the singlet and triplet inter- 
actions, and 



G~ 
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The Q matrix considered here has the following struc- 
ture: Q — Q"f(k — k';ui,ui'), where the pairs ij and af3 
denote replica and spin indices, respectively, while k, k' 
are momenta and to, to' are frequencies. The trace opera- 
tion runs over all these indices. For our purposes, it will 
be sufficient to consider only the zero-temperature limit. 
Let us expand the action for Q = Qo + SQ: 
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where Gq 1 follows from Eq. (||) after replacing Q by Qq. 
The reason for keeping up to order SQ 3 when searching 
for the saddle point will become transparent below. Af- 
ter this expansion, the usual next step is to choose Qo 
such that the term linear in 5Q vanishes; this leads to 
the saddle point used by Finkelstcin. Equivalently, this 
condition on the linear terms can be cast as 
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where Fq 1 ^ is the tree level (or zero-loop) one- vertex func- 
tion. By writing the usual saddle-point equation as a 
condition on the one-vertex Tq 1 ^ (or linear in SQ) poten- 
tial, we are approximating the whole vertex function or 
effective potential T[Q] = V[Q] by the tree level poten- 
tial r [Q] = S[Q]. This leads to the tree level saddle 



point Q "f(q;u),u)') = sgn(o;) £ ?)0 <W S a p Fluctu- 
ations are then parametrized by slow unitary rotations 
Q = C/ _1 (r) Q U(r). The sigma model constrained to 
this manifold has an instability in the triplet interaction 
channel, which signals that the wrong saddle point has 
been chosen. 

We will proceed by considering not the tree level poten- 
tial or action in the derivation of saddle point, but the ef- 
fective potential generated by the fluctuations. There is a 
symmetry breaking not present at the tree (or "classical" ) 
level, but that surfaces at higher loop order ( "quantum" ) 
in the potential, an effect know as the Coleman- Weinberg 
mechanism With this in mind, the correct saddle for 
the interacting disorder electron problem should be found 
at one-loop order: 
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This is the reason why we kept up to the cubic term in 
Eq. (||) . To one-loop order, the cubic term contributes to 
, as shown in Fig|l|. Notice that, if it were not for the 
interactions, this contribution would vanish because it is 
proportional to N — > in the zero-replica limit. The in- 
teraction keeps the contribution alive because it requires 
all replicas in the loop to be the same. 




FIG. 1. One- loop contribution to the one- vertex term in 
the effective potential. 

Contracting two of the SQ's from the cubic term, we 
obtain (spin and replica indices are not shown): 



SV = ^ X! G o(^i' G o(u 2l k 2 ) G (uj 3 , k 3 ) {SQ UlU2 {ki - k 2 ) SQ U j 2UJ3 {k 2 - k 3 )) SQ U3U11 (k 3 - h) (9) 
or, alternatively, 

r i 1) = ^zC [Go(wi,*i)] a G (tu2,k 2 ){SQ^ 2 (k 1 -k 2 )SQ^ 1 (k 2 -k 1 )). (10) 
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The propagator of density fluctuations V = (SQ SQ) fol- 
lows from the quadratic part of Eq. (^|) and depends on 
the saddle Qo through Go- Redefining the frequencies 
and momenta in the sums, one has 



r i 1} =^EE [Go^,fc)] 2 G (c + c',fc + g )P( w ,g) 

uj' k.q 

(11) 

Next, we use the above Ij and find the new saddle solu- 
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tion of Eq. (||). The solution is a matrix Q homogeneous 
in space and diagonal in the replica, spin, and frequency 
indices: Q ?f (q;w,u') = Qq (lj) S q , S uu > 8 aj3 8 i3 ■, where 

we allow for Q (ui) ^ Qq(oj). Let us also assume that 
the electrons interact through a screened (short-ranged) 
interaction. The new saddle-point equation becomes 
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The first term on the right-hand side of Eq. (|12j) is the 
only one present in the non-interacting case. The sec- 
ond term is the one-loop contribution while the third 
one corresponds to the Hartree-Fock approximation. In 
the non-interacting case the solution to the saddle-point 
equations is simply Qq(u>) — sgn(w). For the interacting 
case, we separate the singlet and triplet contributions, 
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where V s and T t are the singlet and triplet coupling con- 
stants, respectively, and a r are Pauli matrices. While the 
singlet channel leads to a shift in the energy band bot- 
tom, a repulsive interaction in the triplet channel causes 
a magnetic instability related to a net spin polarization, 
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similar (but not identical) to the Stoner instability of 
clean, itinerant electrons in the presence of a ferromag- 
netic interaction. In order to capture this effect, we make 
the Ansatz 



Q£H =Csgn(w) -2irA a , 



(14) 



where £ and A a have to be determined self-consistently. 
For simplicity, we set T s = — T t = T/2. 

The contribution from the triplet channel to the T> = 
(SQ SQ) propagator can be written as 



V{u, q ) = - V -±^Df{u, q )Df{^ q ) 
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Here D is the diffusion constant, A a p = A a — A/3, zq = z 
(the frequency renormalization factor), and Z2 = z + vpT . 
It becomes apparent that £ renormalizes amplitudes, sim- 
ilarly to the wave-function renormalization which ap- 
pears in the RG treatment of the problem |],[j],|lO| . By 
dimensional analysis, we expect the bare interaction cou- 
pling constant Y to be replaced by r/£ 2 - we have already 
taken that into account when writing Eq. (|l5|). Thus, the 
contribution to with opposite spins (a (3) in the 
loop is 
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where g = vpD is the dimensionless conductance in 2-D 
and z = [z + Z2)/2. In Eq. ( fL7j ) we have constrained 
M,|A Q/3 | « 1/r. 

With the one-loop contribution Eq. ([l7]), one can now 
replace it in Eq. ( |l2| ) and solve self-consistently for the 
magnetization bandwidth A = Af — Aj_ . After some al- 
gebraic manipulations, we find that the A in the saddle- 
point solution must satisfy 

z [ (2irz) 2 q ( 1 \1 

The upper bound in frequency of the diffusion propagator 
sets the prefactor of the exponential in Eq. (|l8|) . As a re- 
sult, A is proportional to the elastic scattering rate, 1/r, 
rather than the total bandwidth or the Fermi energy. 



Equation ( |l8| ) points to the existence of a non-zero spin 
polarization A > for any positive value of T, provided 
that the dimensionless conductance is finite. Only for 
infinite g we recover the usual Stoner instability char- 
acteristic of clean systems, namely, a ferromagnetic in- 
stability at r > \ jvp. In the particular case of a finite 
system, the same tendency towards spin polarization was 
found by Andreev and Kamenev |ll[ . This suggests that 
the ferromagnetic instability may be a robust property 
of 2D disordered interacting electrons, since their start- 
ing point was rather different than ours. For the finite 
system case, they used an exact basis representation for 
the non-interacting problem, combined with disorder av- 
eraged Hartree-Fock matrix elements. These matrix ele- 
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merits were dressed by the diffusive dynamics to lowest 
order in 1/g and provided a contribution to the magneti- 
zation similar, but not identical, to our one-loop calcula- 
tion with the disorder average carried out from the very 
beginning. 

It has been suggested that the divergence of the triplet 
coupling is connected to the formation of local magnetic 
moments in the system Jl^,|l3|]. Notice, however, that 
the polarization of the saddle-point solution in Eq. JlS| ) 
connects continuously with the Stoner instability in the 
limit of a clean system. This is suggestive of a dirty ferro- 
magnetic state |i~4| where there exists a residual coupling 
between the local moments, causing a tendency towards 
a global ferromagnetic order at T = 0. If the localized 
moments are coupled through, for example, an RKKY 
interaction that has alternating signs, one would expect 
not a ferromagnetic state, but rather something similar to 
a spin glass state [[[5) . There can be two solutions to this 
puzzle. One is that the characteristic size of the local mo- 
ment diverges, and hence the system is in a ferromagnetic 
state (at T = 0) according to the saddle-point solution 
in this paper. The other alternative is that there may be 
another saddle, which in contrast to ours, breaks replica 
symmetry (RS). Let us remind that in spin glass mod- 
els the replica symmetric solution is unstable in the glass 
phase, which is characterized by a non-zero Edwards- 
Anderson order parameter without magnetic order. Our 
saddle-point solution, on the other hand, finds directly 
a magnetic order parameter A ^ 0. In this case, it ap- 
pears unlikely that RS breaking may resolve the issue. In 
addition, any RS breaking solution must be such that it 
recovers the usual Stoner instability in the clean limit. 

Another important issue to consider is the divergence 
of the frequency rescaling factor z in the RG equations for 
the sigma model expanded around the zero-loop saddle. 
In the introduction, we claimed that the divergence of the 
interactions happen at a finite length scale; however, the 
divergence of z would mean that the associated energy 
scale goes to zero. The latter conclusion is not correct 
for the following reason. The coupling constant with the 
leading divergence is the triplet, and when this coupling 
leaves the perturbative RG limit, one can no longer trust 
the flow. The saddle with the non-trivial magnetization 
resolves the leading triplet divergence, and that should 
stop the z divergence as well. 

We would like to conclude the paper by discussing the 
issue of localization in disordered and interacting sys- 
tems. The nature of the triplet channel instability is 
clearly captured by a disordered averaged effective the- 
ory through a non-trivial saddle-point solution with a 
non-zero magnetization. Moreover, the saddle point pro- 
vides a natural starting point of an new RG treatment 
for the localization problem. Once A ^ is taken, it ap- 



pears that the renormalization of the diffusion constant 
D should proceed similarly to the non-interacting case. 
Although a renormalization program should be carried 
out to confirm whether D scales to zero, it appears un- 
likely that starting from the new saddle will lead to de- 
localization in 2-D for weak interactions. This, however, 
would not preclude a correlated (super) conducting state, 
for example, as proposed by Phillips and coworkers S. 
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